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Abstract. The local deformation of two-dimensional Lennard-Jones glasses 
under imposed shear strain is studied via computer simulations. Both the mean 
squared displacement and mean squared strain rise linearly with the length of 
the strain interval A7 over which they are measured. However, the increase in 
displacement does not represent single-particle diffusion. There are long-range 
spatial correlations in displacement associated with slip lines with an amplitude 
of order the particle size. Strong dependence on system size is also observed. The 
probability distributions of displacement and strain are very different. For small 
A7 the distribution of displacement has a plateau followed by an exponential 
tail. The distribution becomes Gaussian as A7 increases to about .03. The strain 
distributions consist of sharp central peaks associated with elastic regions, and 
long exponential tails associated with plastic regions. The latter persist to the 
largest A7 studied. 
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1. Introduction 

The mechanical properties of amorphous materials are important in many structural 
applications and the microscopic mechanisms underlying plastic deformation have 
been the subject of great theoretical interest [H12113] • While deformation of crystalline 
materials can be understood through the generation and propagation of dislocations, 
no analogous defect has been identified in disordered materials. This has a direct 
impact on their mechanical properties since dislocation motion reduces the yield stress 
of crystals by orders of magnitude. In the last decade, it has also been realized that 
amorphous materials share many features with colloids, granular media and other 
systems where motion is jammed or arrested by packing constraints Common 
features are observed in the onset of motion under external forces in such systems and 
understanding their behavior represents a fundamental problem in nonequilibrium 
statistical mechanics. 

Recent simulations and experiments have followed the motion of individual atoms 
or particles to understand how disordered solids deform under shear [21 El O [HI [HI 

MMMMM\MMMM\MMMMMMMMMMMM- One of the 

most studied quantities has been the van Hove function P(Ar; At), which describes the 
probability distribution of the displacements Ar of atoms as a function of time interval 
At. Here Ar is defined relative to a uniformly (affinely) sheared system. In equilibrium 
systems, P(Ar) has been used to determine the time scale on which particles escape 
from local cages [51 [HI 02] and its long time behavior gives the diffusion constant. 
Studies of sheared systems have examined the changes in P(Ar) to determine the 
effect of strain rate on caging times and diffusion [14[ \TE\ |22"1 124] [32] . These studies 
have generally found (Ar 2 ) cx At as expected for single-particle diffusion. In some 
cases the associated nonequilibrium diffusion constant obeys a generalized fluctuation 
dissipation theorem with an effective temperature [TH Q~51 [21] . 

In this work, we use molecular dynamics (MD) simulations to examine the spatio- 
temporal correlations in Ar in athermal systems that are driven by steady shear. 
This models the low temperature behavior of amorphous atomic systems as well as 
disordered granular or colloidal systems where thermal fluctuations are small. As in 
most studies of similar systems [221 [211 [S2] i we find that the long time evolution of 
P(Ar; At) is consistent with the diffusion equation. However, the displacement shows 
strong spatial correlations that span the entire system. Thus P does not describe 
single-particle diffusion and the diffusion constant has a strong dependence on system 
size [32] . The spatial correlations in Ar are associated with plastic deformation along 
slip lines that grow to span the system. The slip accommodated in these plastic zones 
is one or two molecular diameters. This sets a characteristic strain interval between 
their formation that decreases with system size. Spatio temporal correlations in the 
location of the slip lines are studied through images and the evolution in the functional 
form of P with strain. 

The long-range spatial correlations in Ar reflect the fact that large regions deform 
elastically. To measure the local magnitude of non-affine deformations, we evaluate 
the vorticity in the displacement field to. Other measures of the relative particle 
displacement have been considered in previous work, but lu has the advantage that it 
allows both the magnitude and direction of the displacements to be quantified. This 
allows strong directional, spatial and temporal correlations in strain to be identified. 

The distributions of Ar and u> provide additional information about correlations. 
Over short intervals, both distributions show exponential tails covering up to eight 
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orders of magnitude. The strain characterizing the exponential decay for us is close 
to the the macroscopic yield strain of the material. As the interval of time or 
strain increases, P(Ar) becomes more isotropic and converges to a Gaussian at strain 
intervals of order 3%. The distribution of us retains an exponential tail out to much 
larger strains, indicating that particles are still not undergoing single-particle diffusion. 

2. Protocol and definitions 

We perform classical 2D molecular dynamics simulations of a bi-disperse mixture that 
is chosen to inhibit crystallization 29J. There are Nl large (L) particles and N$ 
small (S) particles with N^/Ns = (1 + \/5)/4. The interaction between particles of 
type I and J is a Lennard Jones (LJ) potential: Uu — 4e [(r/au)~ 6 — (r /au)~ 12 ] 
where e sets the characteristic energy scale and 07 j reflects the size of / and J. We 
set gll = 1-Octo, ass = 0.6<7oi ancl o~ls = 0.5(oll + ass) = 0.8<7o an d measure all 
lengths in units of Co- To speed calculations, the potential is truncated smoothly using 
a fourth order polynomial that matches the LJ potential and force at ri n = 1.2<t/j and 
takes both to zero at r ou t = l-5crjj. All particles have mass m and the characteristic 
time is t = yj ma^/e. 

The equations of motion are integrated using the parallel MD code LAMMPS 
|33j with time step St = 0.0056t. The system temperature is maintained at zero by 
applying a locally Galilean- invariant Kelvin damping mechanism [34 . Neighboring 
particles exert equal and opposite drag forces on each other which are proportional 
to their relative velocity. The value of the damping constant was set to QAbm/r. 
With this choice, modes with wavelengths at the particle scale are nearly critically 
damped and damping decreases at longer lengths. The limit of zero temperature is of 
fundamental interest in atomic systems and thermal motion can generally be ignored 
in systems where the particles represent larger entities such as grains or colloids. 

An amorphous state of the system is prepared following the protocol in previous 
studies [30] . Periodic boundary conditions with periods L x and L y are used. A random 
configuration with a particle area fraction of 0.85 is equilibrated for 20, 000 time steps 
using a soft repulsive potential to ensure particles do not overlap. The interaction is 
then switched to the LJ potential, and the damped system is run for 10, 000 timesteps 
with fixed area. The system is next allowed to expand toward a zero pressure state 
under an isotropic Nose-Hoover barostat [33] with a time constant of 220r and a 
damping parameter of 45 r _1 for 10,000 timesteps. A final quench is run for 10,000 
timesteps with a fixed cell size prior to shearing. 

Equilibrated systems are compressed along the y direction at constant true 
strain rate 7 and expanded along x to maintain constant area: L x = L X Qe +Jt , 
L y = Lyoe - ^*. Unless noted, the results below are for about 1.6 x 10 6 particles 
with L x0 — L y0 w lOOOoo: and 7 = 0.45 x 10~ 5 r -1 . Studies showed that this strain 
rate was slow enough to be in a quasistatic regime where results are rate-independent. 
For this reason, results will be expressed in terms of the interval of strain A7 over 
which they were obtained, rather than in terms of time interval. For small strains 
the rate of plastic deformation increases with strain and may be sensitive to initial 
conditions [HI [201 ES] . For strains greater than 6% the behavior becomes independent 
of strain and initial state. Data presented below are all from this steady state regime. 

The non-affine component v of a particle's velocity v is the deviation from the 
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velocity associated with a uniform strain: 

v x — v x - jx(t) , . 

Vy = Vy + 7?/ (r) 

Integrating over time gives the net non-affine displacement Ar = (Ax, Ay) over any 
finite strain interval. Note that this is not the same as taking the difference between 
the final position of a particle and the position corresponding to an affine displacement. 
The reason is that a small non-affine displacement at the beginning of a time interval 
is amplified by subsequent affine flow and would have the same impact as a much 
larger non-affine motion at the end of the interval. A similar effect produces an extra 
Taylor diffusion in sheared fluids [36j . From now on we shall refer to the non-affine 
displacement simply as "displacement" and use "total displacement" when the affine 
component is not subtracted. Note that the non-affine displacement averaged over the 
periodic cell is zero. This follows because the equations of motion in the deforming cell 
|33j conserve the sum over particles of the peculiar momentum rav and all particles 
have the same mass. 

The results presented below show that the displacement field has long-range 
spatial correlations. A more local measure of deformation is provided by the spatial 
derivatives of Ar, which are obtained in the following way. A Delaunnay triangulation 
is found for an initial reference configuration, with vertices at the center of all particles. 
The derivatives diArj on each triangle are then calculated from the finite difference 
between displacement vectors at the three vertices. To make plots and evaluate 
probability distributions, the derivatives are sampled on a fine, rectilinear grid with 
spacing ~ cr /3. An important property of this piece-wise constant definition of diArj 
is that its integral along an arbitrary path connecting any two particles gives the exact 
difference between the displacement vectors of the two particles. This guarantees that 
mean strains on large scales are equal to the spatial average of local values. 

Because the area is fixed and the system is densely packed, changes in density 
are small. As a result, we find that the local divergence of the displacement field 
is negligible. Most of the deformation takes the form of local shear strains that are 
quantified by the curl of the displacement field or vorticity: lo = d y Ax — d x Ay, 
which represents the antisymmetric part of the derivative tensor. For two dimensional 
systems the curl is a vector perpendicular to the plane and thus invariant under 
rotations within the plane. Note that it is conventional to associate the magnitude of 
shear with a symmetric scalar invariant, }^\J (d x Ax — d y Ay) 2 + (d x Ay + d y Ax) 2 , the 
so-called deviatoric strain. Such symmetric definitions of the strain are completely 
blind to the vorticity of the displacement which we show below to be quite important 
in our system. In particular we will see that the sign of u> is strongly correlated 
with the orientation of lines along which shear occurs. As the magnitude of the 
strain increases, nonlinear terms in the definition may become important. The 
vorticity can be generalized to a rotation tensor R^ defined by the decomposition: 
Sij + diArj — SikRkj where Sy is a symmetric tensor which describes pure axial 
stretching. In the limit of infinitesimal strains, u> and the rotation angle associated 
with Rij become identical. For finite strains, we find qualitatively similar results for 
the distribution of u and the rotation angle. 




Figure 1. Grayscale plots of the projections of a typical non-affine displacement 
field u onto different directions. Top left: Ax. Top right: Ay. Bottom left: 
(Ax + Ay)/%/2. Bottom right: (Ax - Ay)/y/2. Here A7 = .002, the x-axis is 
horizontal, and the gray scale is linear, ranging from — ctq (black) to +cro (white). 

3. Spatial Variation of Displacement and Strain 

In figure [TJ we plot projections of a typical displacement field for A7 = 0.002 
onto horizontal, vertical and diagonal directions: Air, Ay, (Ax + Ay)/\/2, and 
(Ax— Ay) / -s/2. All projections show strong spatial correlations over scales comparable 
to the system size, but different projections highlight different features. For example, 
the plot for the component proportional to Aa; — Ay is dominated by a region of 
strong gradient extended roughly along y ~ ~x while the Air + Ay component has a 
few smaller features extended along y — x and centered roughly about the large feature 
in Aa; — Ay. The displacements Aa; and Ay show a superposition of all these features. 
Note that some of the features extend around the periodic boundary conditions. 

The strong spatial correlations evident in Fig. Q] show that most of the 
displacements of atoms that contribute to P(Ar; A7) are very far from independent 
random walks by each particle. The correlations result from shear along spatially 
extended slip lines. The slip lines show up as lines of rapid gradients in Fig. [TJ Because 
they reflect a shear displacement, the amplitude of the displacements is largest for the 
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Figure 2. An idealized schematic representation of the (top) total and (bottom) 
non-affine displacements parallel to a slip line as a function of distance away 
from the line. The change a in the total displacement across the periodic cell 
is accommodated in a plastic zone of width h where the displacement changes 
rapidly, u> PS a/h. The total displacement in the surrounding clastic regions is 
nearly constant. Subtracting the affine displacement associated with a strain A7 
leads to a small gradient, u) ~ A7, in the elastic regions of plots of the non-affine 
displacement. 

projection along the slip line. Thus the slip line that spans the entire system along 
x — y shows up most clearly in Ax — Ay, while a number of smaller slip lines along x + y 
are emphasized in Aa; + Ay. The reason that slip lines are at an angle of roughly 45° 
from the compressive axis is that the stress tensor has the largest shear (off-diagonal) 
components for these orientations. Q Put another way, shear along these lines is most 
efficient in simultaneously contracting the system along the compressive direction and 
expanding it in the tensile direction. 

Figure [2] illustrates (top) the total displacement and (bottom) the non-affine 
displacement along a slip line as a function of the distance away from the slip line. The 
change a in the total displacement across the periodic cell is accommodated almost 
completely in a plastic zone of width h. The total displacement in the surrounding 
elastic regions is nearly constant. Subtracting the affine displacement associated with 
a strain A7 leads to a small gradient u> « A7 = a/L in the elastic regions of plots 
of the non-affine displacement. For h « L the gradient in the plastic zone is much 
higher uj ss a/h. 

Traces through Fig. [Tj and similar plots are roughly consistent with Fig. [5] In all 
cases there are large regions where the total displacement is nearly constant, leading 
to a gradual slope uj ~ a/L in the non-affine displacement. The amount of strain 
accommodated in a plastic zone varies somewhat with the strain interval, particularly 
for slip lines that do not span the system and where there are multiple slip lines in a 
given strain interval (Fig. Q] bottom left). The edge of the plastic zone is also more 
diffuse than in Fig. [21 leading to a rounding of the peaks in the non-affine displacement. 
This rounding varies along the slip line as is evident from the fluctuations along the 
primary ridge in the bottom right panel of Fig. [TJ These variations have important 

f Note that the periodic boundary conditions may favor slip line orientations along low order periods 
such as the diagonals, but we find the angle remains near 45° to the compressive axis even when the 
cell has deformed significantly. 
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implications for the distributions discussed in the next section. 

The form of plots like Fig. [1] depends on strain, strain interval and system size. 
At very small strains, the entire system responds elastically and the displacements are 
small. As the strain increases, the stress and rate of plasticity rise. For strains of 6% 
and greater we find that the system is in a nearly steady state. In this regime, plastic 
deformation tends to occur in a series of rapid "avalanches" with a wide range of 
sizes that will be discussed further in future work [37]. Between these avalanches the 
displacement field may be small, but the strain interval between avalanches decreases 
rapidly with system size. For the system size considered here, plastic zones that span 
a large fraction of the system are always observed for A7 > 0.002 and this is the 
reason we focus on intervals of this order. Many avalanches occur during this strain 
interval and they cluster spatially to produce the slip lines evident in Fig. Q] 

The evolution of the spatial organization of plastically deformed regions is 
illustrated in Fig. [3] The left hand side shows the magnitude of the non-affinc 
displacement, Ar, rather than individual components. It is these magnitudes that 
are used to construct the widely studied P(Ar;A7). The right hand side shows 
to, which quantifies the amount and direction of local shear. The successive panels 
moving downwards in Fig. [3jshow the changes over progressively longer strain windows 
starting from the same initial configuration as in Fig. [TJ 

The families of slip systems along both diagonals contribute equally to Ar and are 
visible in the left hand panels. However, taking the magnitude of the displacement 
obscures the nature of the atomic rearrangements. The slip lines that appeared as 
sharp steps in Ax + Ay and Ax — Ay appear now as regions of low Ar surrounded 
by halos of higher Ar. The origin of this behavior is evident from Fig. [2] The peak 
displacements occur at the outer edges of the plastic zone at a distance of order h/2 
from the center of the slip line. The magnitude drops rapidly to zero at the slip line 
and much more slowly in the elastic region outside. 

A better measure of the local plastic response is provided by uj. The oj field 
is clearly localized on the regions of strong shear and much smaller away from these 
regions. From Fig. [2] we expect the shear strain to be of order A7 in the elastic regions 
and much higher in the plastic regions that accommodate almost all the deformation 
of the cell. The numerical values of \oj\ are consistent with this expectation. Note that 
the sign of u> carries important information about the direction of slip. In particular, 
the two families of slip planes produce opposite signs, with a clockwise rotation along 
lines running parallel to x — y and a counterclockwise rotation on perpendicular slip 
lines. Thus plots of ui highlight the correlations along individual slip lines. 

The successive panels moving down Figure [3] show how displacements add as 
the strain interval increases. Over short intervals A7 < 0.002 we find that the 
displacements tend to organize along system-spanning lines. This allows the entire 
displacement to be accommodated in plastic zones with minimal energy stored in 
deformation of elastic regions. Subsequent slip lines appear to occur at random 
locations rather than continuing to nucleate along the same path to produce a 
persistent shear band (e.g. as in reference [38]). One way of quantifying this is 
through the autocorrelation function for displacements over successive intervals. We 
find that this correlation function drops rapidly. By A7 = 0.003 the correlations have 
dropped by about 2 orders of magnitude and are comparable to the noise. This is 
consistent with Fig. [3] where we see a growing number of independent slip lines for 
A7 > 0.004. A superposition of randomly located slip lines is consistent with the 
diffusive growth in Ar and oj discussed in the next section. 
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The strain interval over which deformations are correlated depends on system 
size and seems to be related to a characteristic amplitude of the slip across plastic 
zones. Examination of many plots like Fig. [1] shows that a varies from about the 
mean particle diameter 0.8cro in the central region of shorter slip lines to about twice 
this value in a system spanning slip line like that in the bottom right of Fig. [TJ It is 
perhaps not surprising that once a is of order a particle size it is energetically favorable 
to have plastic deformation everywhere along the slip line rather elastic deformation. 
Note that the strain interval needed to produce a displacement a is of order a/L. 
This implies that throughgoing faults should occur over A7 ~ 0.001 to 0.002 in our 
systems. Our observations are consistent with this and analysis of the avalanche size 
distribution also shows that the largest events are spaced by A7 ~ 0.002. 

The plastic regions in Fig. [T] and on the left of Fig. [3] have a significant 
width h ~ 50er. Thus the mean shear strain across these regions is only of order 
a/h ~ 0.02 — 0.04 for A7 = 0.002. While this is an order of magnitude larger than the 
mean strain A7, the peak strains on the right of Fig. [3] are higher still. Each broad 
slip zone contains a number of much sharper features with strains that are an order 
of magnitude higher (\oj\ = 0.2 -0.4). The width of these smaller features is only a 
few particle diameters. More detailed studies of the time dependence show that the 
slip zones are formed by a large number of avalanches that are spatially correlated to 
produce the system spanning lines in the figures. The long-range spatial correlations 
in the strain field are analyzed in Ref. [39) and we will focus below on the temporal 
evolution of the probability distribution of displacement magnitudes and u). 

Comparing results for different system sizes indicates that h increases roughly 
linearly with Lq while the width of the high strain regions remains of order the particle 
diameter. It is interesting to note that the width of fault regions in earthquake systems 
also tends to scale in rough proportion to the length [40] . Studies of different system 
sizes also indicate that the slip amplitude along throughgoing faults is always of order 
one to two particle diameters. If a is independent of system size, then the strain 
interval needed to produce a slip line should be inversely proportional to system size. 
There should also be a slower decay with A7 in the displacement autocorrelation 
function. Results for smaller systems (Lq — 250 and 500) are consistent with these 
predictions. 

4. Probability distributions of displacement and strain 

Figure 2] shows the probability distribution for finding a given displacement over 
successive strain intervals with A7 = 0.002. The distribution P(Ar) is normalized 
so that the integral over all vectors Ar is unity. With this definition the distribution 
goes to a constant as Ar — > 0. Note that the distribution varies significantly from 
interval to interval. This reflects the large scale structure of the slip lines discussed 
above. In contrast, the distribution of the more local variable oj is nearly the same for 
all intervals. Its average behavior is discussed below. 

For an ideal system-spanning slip line with a constant magnitude a and h « L, 
traces like Fig. [2] (bottom) would consist of two straight segments of constant slope 
A7 and a/h, respectively. The distribution P(Ar) would be constant from to a/2 
and zero for larger displacements. Given the range of a discussed above, one would 
expect plateaus extending to Ar ~ 0.4<7o to 0.8<7o in Fig. 2J The curves show a 
crossover from gradual to rapid decrease at Ar in this range, but do not show true 
plateaus. Analysis of displacement fields like those in Figs. Q] and [3] show that a varies 



Figure 3. Left: The magnitude, Ar, of the (non-affine) displacement fields. 
Gray scale is linear ranging from (black) to +2.5<ro (white). Right: The curl, 
lu, of the (non-affine) displacement fields. Gray scale is linear ranging from — .25 
(black) to +.25 (white). Displacements arc computed relative to the same initial 
state as in Fig. [T]and over growing strain intervals A7 = .002 (top), .004 (center), 
and .008 (bottom). 



Evolution of displacements and strains in sheared amorphous solids 



10 




At 

Figure 4. P(Ar) for several consecutive intervals with A7 = .002. The average 
over all intervals from 7 = 12% to 24% is shown in heavy black. 

along the slip line. As discussed above, there is also smearing along the edge of the 
plastic region that rounds the peaks in traces like Fig. [5] Interactions between the 
two families of slip planes also alter the distribution. Together theses effects lead to 
a reduction in the range of Ar where the distribution for each interval is constant, 
and an exponential tail in the probability distribution at large magnitudes. Averaging 
over all intervals between 12% and 24% produces the thick solid line. The crossover 
to the exponential tail is broadened by averaging over intervals with different values 
of a and the exponential tail becomes more pronounced. 

Figure [5] shows how the mean squared displacement and strain increase with A7. 
All results are averaged over strains between 12% and 24%. As would be expected 
for a diffusive process, both (Ar 2 ) and (oj 2 ) rise linearly with A7. The fit lines 
have a small negative offset at A7 = 0. This reflects the correlation in the increase 
in both quantities at small strains, which cause them to rise quadratically with A7 
at small strain intervals (A7 < 0.001). While the linear rise of (Ar 2 ) at large A7 is 
consistent with diffusion, this only implies that the displacements of any given particle 
in successive strain intervals are uncorrelated. As shown above, there are strong spatial 
correlations that imply one should not view the slope of Fig. [5jleft) as a single particle 
diffusion constant. A detailed examination of the probability distributions also reveals 
some subtle longer- lived correlations in displacement. 

Figure [6] shows the probability distribution of displacements averaged over initial 
states and scaled by \AOOA7 to remove the mean increase in magnitude shown 
in Fig. [5] For the smallest interval, the distribution has an exponential tail that 
extends over eight orders of magnitude. It is interesting to note that the displacement 
distributions following individual plastic events or avalanches have also been found to 
be exponential [22l[32] . If successive events were strictly decorrelated, the distributions 
for longer intervals would be given by convoluting the distributions for short intervals 
and would converge to a Gaussian. The distributions do seem to become more 
Gaussian as A7 increases, but there is a systematic decrease in the probability of 
large events for A7 = 0.008. This suggests that the behavior is still correlated, with 
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Figure 5. Variation with A7 of the mean squared (left) displacement (Ar 2 ) and 
(right) strain (a; 2 ). Straight-lines are linear fits through the data points, which 
have statistical errors that are slightly smaller than the symbol size. 



large displacements less likely to occur in the same location as in previous intervals. 
By A7 = 0.032 these correlations have nearly disappeared and the data are consistent 
with a Gaussian distribution. It is interesting to note that at this point the product 
of the number of throughgoing faults (~ 0.032/0.002 = 16) and the fraction of the 
system affected by each (h/L ~ 50/1000) is of order one, implying that slip lines must 
revisit the same regions of space. 

Figure [7] shows the time averaged u) distributions for various A7. The 
distributions consist of central peaks that cross over to exponential tails. The 
exponential tails persist even at the largest A7. Roughly speaking, regions of the 
material with small strains away from any shear zones (the gray regions in the right 
column of Figure [3]) contribute to the "elastic" central peaks, while the strongly 
sheared material in the core of the shear zones (the black and white regions in the 
right column of Figure [3]) contributes to the "plastic" exponential tails. To first order, 
weight is removed from the central peak and redistributed to the tails as the amount 
of plastic strain increases. 

To test this picture quantitatively, the u distributions are scaled by A7 and 
replotted in figure [5J This rescaling comes close to collapsing the data in the tails for 
the smallest A7, but becomes less satisfactory for increasing A7. Direct convolution 
of data from small intervals shows similar deviations from probabilities for larger 
intervals. This implies that there is some additional correlation over the range of 
A7 considered here. Fits to an exponential decay, P(uj) ~ exp(w/w*), show a slight 
increase in to* with A7 (inset). This increase implies that increasing A7 does not 
merely increase the size of the plastic region and thus the total weight in the tails. In 
addition, there is a slight increase in the relative weight of large magnitude regions 
that must reflect a correlated increase in strain in regions of large strain. 

It is interesting to note that the value of w* has a natural physical interpretation. 
The vorticity w = d y Ax — d x A y while the shear strain is 0.5 * [d y Ax + d x A y ]. Thus 
uj* =0.1 corresponds to a shear strain of about 5%, which is close to the yield strain 
|19j As noted above, the regions of large u> have strains much bigger than the mean 




Figure 6. Average displacement distributions, P(Ar;A , y), for A7 = 0.001 
(dot-dashed), 0.002 (dotted), 0.008 (dashed), and 0.032 (solid). To compensate 
for the increase in width with A7, the curves are plotted as a function of 
s = Ar/ao (IOOA7) ' 5 . The thick (red) dotted line is a one-parameter fit of 
the A7 = .032 data to a Gaussian. Averages are over all intervals from 7 = 12% 
to 7 = 24%. 




Figure 7. The time averaged u) distributions, P(lo; A7) for A7 = 
.001,.002, .004, .008, .016. 
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Figure 8. Distribution of w scaled by A7 in percent for A7 = 
.001, .002, .004, .008, .016. Curves increase with A7 at large u. Straight lines 
show exponential fits, P ac e _ "/"« j and the inset shows w* as a function of A7. 

strain in the plastic region (a/h) but small enough that there are several in the plastic 
zone that combine to produce the total displacement. The probability distribution 
suggests that the yield strain sets the characteristic scale for these high strain regions. 

Another important parameter is the strain interval required to produce plastic 
deformation throughout the entire system. One can estimate this from the integral 
under the exponential tails in P(u>) at small strains, and the assumption that the 
weight grows approximately linearly with A7 as in Figure The result is that the 
regions of plastic strain should cover the system after a strain of about 10%. This is 
larger than the estimate of 3% strain required for plastic zones to cover the system, 
but as noted above, only a fraction of the plastic zones has the large strains indicative 
of local plasticity. 

5. Discussion 

Although ours is, to our knowledge, the first study of strain distributions in simulations 
of sheared systems, several studies have looked at displacement distributions in 
athermal systems. All considered simple shear, v x = 77/, rather than the pure 
shear strain studied here. Some simulations used quasi-static simulations with 
small strain steps [HJ [2H 132], while others used a constant slow shear rate as 
in our simulations [21 [TS] SI]. The quasi-static studies all employed Lennard- 
Jones potentials and imposed simple shear with Lees-Edwards periodic boundary 
conditions [211 [32] or rigid walls parallel to the flow direction [22]. Ono and co- 
workers [TU [15] simulated steady shear of wet foams using Durian's bubble model [42] , 
while Radjai and Roux simulated granular media with dissipation entering through 
friction at particle contacts [41 j . Only the latter work and our simulations include 
inertia. Despite the wide range of protocols and interactions, almost all studies find 
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diffusive behavior for the non-affine displacement. |j The only exception is the work of 
Radjai and Roux [UJ who found (Ay 2 ) cx A7 18 . It seems very important to determine 
the cause of this striking difference. The only unique feature of these simulations 
appears to be the incorporation of friction forces, which are not generally thought to 
lead to such qualitative changes in behavior. 

The value of the effective diffusion constant D e ff = (Ar 2 )/A r y can be estimated 
from our observations of slip lines and the simple model in Fig. [2] This displacement 
field gives a constant distribution of displacements from —a/2 to a/2, implying a mean- 
squared displacement of a 2 /12. If system-spanning slip lines produce independent 
displacements at strain intervals of a/L, then (Ar 2 ) w (A7 /(a/ L))a 2 / 12 or -D e // = 
La/12. Using the value of D e ff — 57cr§ from Fig. [5j one finds a s» 0.7<7rj. Improving 
the estimate of the mean-squared displacement from each slip line leads to larger values 
of a since the actual distribution of displacements (Fig. [4]) is cut off before a/2. We 
conclude that the measured diffusion is consistent with the direct observation that slip 
lines have amplitudes of one to two particle diameters and contribute independently 
to the displacement. 

This model for the effective diffusion constant implies a linear scaling with system 
size. We have not yet examined this scaling directly, but recent studies of a similar 
model by Lemaitre and Caroli |32j showed a pronounced increase in D e f f with system 
size. Data for the two largest systems, L = 20 and 40, appear consistent with linear 
scaling with L. Results for L = 10 are larger than expected, but there may be 
significant finite size effects in these systems. Note that we find h/L ~ 1/20 so system 
sizes with L < 20 could not have the plastic zone occupying the same fraction of the 
system. Extrapolating Lemaitre and Caroli's result [jf, for L = 40, D e ff — (2.2 ± .1)<J 2 
to L = 1000 yields D e ff w (55 ± 2.5)cr 2 , which is consistent to our measured value 
of 57fig. The magnitude of A7 needed to observe diffusive behavior is also roughly 
a factor of 25 smaller in our systems. These results provide a dramatic check of the 
linear scaling with system size in diffusion, although we caution that these comparisons 
should only be viewed qualitatively, as the impact of different shearing geometries has 
not been taken into account. 

Lemaitre and Caroli discuss the system-size dependence of their results in terms 
of displacements produced by independent avalanches rather than independent slip 
lines. Earlier studies of individual plastic events showed that the rate of large 
avalanches scales with the system length [27l [29]. Based on the assumption that 
the total strain A7 was accommodated by these avalanches, Maloney and Lemaitre 
concluded that the displacement a a associated with each avalanche was independent 
of system size and of order 0.1er . Extrapolating their results [2TJ [29] for L ~ 50 
to our system sizes would give large avalanches spaced by A7 = 0.0002, an order of 
magnitude more frequent than the rate of system-spanning slip lines in our simulations. 
This is consistent with our conclusion that the slip lines in Figs. 1 and 3 are not 
the result of individual avalanches, but the correlated displacement produced by 
multiple avalanches. Assuming that each avalanche contributes independently to the 
displacement would give too small a value of D e ff. 

§ For simple shear this is most easily measured in the displacement perpendicular to the flow 
direction x, (Ay 2 ) <x A7. 

|| Note that in reference |32| . the imposed homogeneous simple shear flow field is defined as v x = yy, 
v v = 0. The associated symmetric strain rate tensor associated with this flow field has a magnitude 
of 7/2, so one needs to scale the values of (Ay 2 ) / A7 from reference [32] by a factor of 2 to make a 
direct comparison to our results for pure shear geometry, v x = 72:, v y = —'yy. 
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System size dependence of -D e // could be crucial in analyses which invoke the 
notion of an effective temperature based on an effective Stokes-Einstein relationship 
between D e f t and the viscosity ij [T51 121] . The product of diffusion and viscosity is 
related to the temperature in equilibrium systems. In systems with a yield stress, the 
viscosity diverges like the inverse strain rate, r\ ~ At/ Aj, at low strain rates where the 
stress reaches a plateau. As we've seen here, one also expects the effective diffusion 
constant to grow linearly with the strain rate: Ar 2 / 'At = (Ar 2 / A r y)(A'y / At). This 
means that the product of the viscosity and diffusivity should become constant at 
small enough strain rate: (Ar 2 / At)r) ~ Ar 2 / A7. However our results imply that this 
constant should have a dramatic size dependence, and one must proceed with caution 
to utilize it as an effective temperature. 

The scaling of -D e // would be very different if the correlations in displacement 
did not span the system. There is some evidence that system spanning slip lines are 
suppressed by increasing temperature jT5J [5U] . They may also be suppressed by the 
surrounding fluid in experiments on colloidal systems [43j . Dissipation through flow in 
the fluid is reduced by maintaining a linear flow profile that may inhibit localization 
of the plastic deformation. Studies of both effects would be of interest. 

The most important message in our results is that one must go beyond the 
second moments of the displacement and strain distributions to get insight into the 
microscopic processes at play in sheared systems. Displacement distributions may be 
nearly Gaussian, scaling with \J At, while at the same time, there may be strong spatial 
correlations and heterogeneities in the system. The strain distributions can remain 
far from Gaussian, even while the displacement distributions become nearly Gaussian. 
Such behavior provides a clear signature of spatial organization. Furthermore, while 
individual particle trajectories may seem jerky with smooth motion interrupted by 
abrupt jumps, these jumps should not necessarily be interpreted as a particle escaping 
the cage formed by its neighbors (as is often assumed to be the case), but could just as 
easily correspond to coherent displacement of the neighborhood with very little relative 
motion of particles. In the future, it will prove interesting to repeat the analysis 
here focusing on strain in addition to the more commonly studied displacement in 
thermal (e.g. super-cooled liquids) or brownian(e.<7. colloidal suspensions) systems 
on approach to the jamming transition. This kind of measurement could serve as an 
important probe of the co-operative nature of the dynamics. 
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